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We extend existing work on the propagation of ultra-high energy cosmic rays in extragalactic 
magnetic fields to a possible component of heavy nuclei, taking into account photodisintegration, 
pion production, and creation of pairs. We focus on the influence of the magnetic field on 
the spectrum and chemical composition of observed ultra-high energy cosmic rays. We apply our 
simulations to the scenarios proposed by Anchordoqui et al., in which Iron nuclei are accelerated in 
nearby starburst galaxies, and show that it is in marginal agreement with the data. We also show 
that it is highly unlikely to detect He nuclei from M87 at the highest energies observed ~ 3 10^° eV 
as required for the scenario of Ahn et al. in which the highest energy cosmic rays originate from 
M87 and are deflected in a Parker spiral Galactic magnetic field. 



I. INTRODUCTION 

The origin of ultra-high energy cosmic rays (UHECR) 
is one of the major open questions in astro-particle 
physics. Data from the Fly's Eye experiment |Q] suggest 
that the chemical composition is dominated by heavy 
nuclei up to the ankle {E ~ 10^®-^ eV) and then pro- 
gressively by protons beyond, while other data |^ may 
suggest a mixed composition of both protons and heav- 
ier nuclei. The fact that present experiments do not give 
a clear answer to the question of chemical composition 
of primary particles motivates to test scenarios with a 
heavy component. 

Nucleons cannot be confined in our galaxy at energies 
above the ankle; together with the absence of a corre- 
lation between their arrival directions and the galactic 
plane, this suggests that if nucleons are primary particles 
they should have an extragalactic origin. At the same 
time, nucleons at energies above ~ 4 x 10^^ eV interact 
with the photons of the Cosmic Microwave Background 
(CMB) by photopion production; this would predict a 
break in the cosmic ray flux, the so-called GZK cut-off §, 
and the sources of UHECR above the GZK cut-off should 
be nearer than about 50 Mpc. The GZK cut-off has not 
been observed by the experiments such as Fly's Eye 0, 
Haverah Park ||, Yakutsk ||, and AGASA However, 
currently there seems to be a disagreement specifically 
between the AGASA ground array which detected 
about 10 events above 10^° eV, as opposed to about 2 ex- 
pected from the GZK cut-off, and the HiRes fluorescence 
detector which seems consistent with a cut-off [ p^ . 
The resolution of this problem may have to await the 
completion of the Pierre Auger project which will 
combine the two existing complementary detection tech- 
niques. 

In the acceleration scenario, UHECR can achieve these 
extreme high energies by acceleration in shocked magne- 
tized plasmas in powerful astrophysical sources, such as 
hot spots of radio galaxies and active galactic nuclei @. 

Attributing sources to the highest energy events is 



complicated by the lack of observed counterparts |T^, ^ . 
A possible explanation is the existence of large scale 
intervening magnetic fields with intensities B ^ 0.1 — 
1 /iC [p^ , which would provide sufficient angular de- 
flection even for high energies and could explain the 
large scale isotropy of arrival directions observed by the 
AGASA experiment Q as due to diffusion. In this frame- 
work, the clusters of events seen by the AGASA and 
Yakutsk experiments |7j are interpreted as due to fo- 
cussing of the highest energy cosmic rays in caustics of 
the extra-galactic magnetic fields, as originally suggested 
in Ref. |^ (see also Ref. for nuclei propagating in the 
Galactic magnetic field and Ref. for recent detailed 
analytical studies). Indeed it has been realized recently 
that magnetic fields as strong as ~ l^G in sheets and 
filaments of large scale structures, such as our Local Su- 
percluster, are compatible with existing upper limits on 
Faraday rotation |lj, |l8[ [l|] . 

Heavy nuclei as UHECR primaries are interesting in 
two ways in this context : they can be accelerated more 
easily to high energies, as the maximal acceleration en- 
ergy a particle can achieve depends linearly on its charge 
Ze, and, in addition, the increased deflection (also pro- 
portional to Ze), could explain more easily the absence 
of correlation between the arrival direction of the events 
and the nearest powerful astrophysical objects. However, 
even in this case there is a limit on the distance to the 
source because of photodisintegration processes due to 
the interaction with infra-red and CMB. 

The study of the propagation of heavy nuclei in the 
absence of magnetic deflection has been treated in some 
detail in the literature. The pioneering work of Puget, 
Stecker and Bredekamp (PSB in the following) which 
included all energy loss mechanisms, has been recently 
updated ^ to take into account new empirical esti- 
mates of the infrared background density of photons Q 
which are about one order of magnitude lower than used 
by PSB. 

In this paper we study the propagation of a distribu- 
tion of heavy nuclei in a stochastic magnetic field, in- 
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eluding all relevant energy loss processes. Our numer- 
ical simulations allow to treat in a consistent way the 
interplay between magnetic deflection and photodisinte- 
gration losses. We also keep track of the propagation of 
all nuclcon secondaries produced in photodisintegration 
events, and propagate these secondaries in the magnetic 
field. These effects had not been considered in previous 
studies of UHE nuclei propagation. In particular, we fo- 
cus here on the influence of the magnetic field on the 
observable UHECR spectrum and its chemical composi- 
tion. In contrast to the sky distribution, these quanti- 
ties are not significantly infiuenced by Galactic magnetic 
fields which we therefore neglect. As will be seen in what 
follows, a relatively strong magnetic field {B > 10^^ G), 
i.e. such that UHECR of low energies diffuse, modify 
by its presence the chemical composition and the energy 
spectrum recorded at a given distance. This is due to the 
effect of diffusion, which increases the local residence time 
differentially with energy, as well as the effective length 
traveled hence the photodisintegration probability. The 
interplay between these effects is rather complex, and the 
output spectrum and chemical composition thus depend 
on several parameters such as the maximum injection en- 
ergy, injection spectral index, linear distance and initial 
chemical composition. Due to the rather high dimension- 
ality of the parameter space, we will show results for fixed 
values of the maximum injection energy i?max — 10^^ eV 
and spectral index dn/dE oc E^^, at the expense of gen- 
erality, and discuss how the conclusions would be modi- 
fied for other values of these parameters. 

The paper is organized as follows: in section II we 
describe the propagation of UHE heavy nuclei, in section 
III we describe our numerical simulation, in section IV 
we present our results, in section V we apply our results 
to test the validity of some recent models, and in section 
VI we conclude. 



number of nucleons emitted. The Lorentz factor of the 
nucleus is given by F = E/{Ampc'^) of the nucleus, e 
and e' are the background photon energy in the observer 
frame and in the rest frame of the nucleus, respectively, 
and n(e) is the photon density of the ambient radiation. 

The range of energies for the photodisintegration pro- 
cess, in terms of the photon energy e' in the rest frame 
of the nucleus, splits into two parts. The first contri- 
bution comes from the low energy range up to 30 MeV, 
in the Giant Dipole Resonance region, where emission 
of one or two nucleons dominates; the second contribu- 
tion comes from energies between 30 MeV and 150 MeV, 
where multi-nuclcon energy losses are involved. Above 
150 MeV, following |20, 2^, we approximate the pho- 
todisintegration rates by zero. This energy corresponds 
to the threshold for photopion production, and we in- 
clude this loss by using the cross-section of nucleon pho- 
topion production scaled by the geometrical factor A^^^ . 
Note that the energy carried away by a pion in such 
an interaction is ~ 20% of the interaction nucleon en- 
ergy, hence only ~ 20%/^ of the primary nucleus. The 
threshold for photopion production is also increased to 
~ 4 10^^ X ^eV, and therefore pion production is only 
important for nuclei up to A ^ 4. 

Returning to photodisintegration, the lower limit of 
the integral on de' , in Eq. (m) was approximated by 2 
MeV for all reaction channels in PSB. We prefer to fol- 
low the approach of Ref. ||2^, with different thresholds 
for emission of one, two and multiple nucleons, for dif- 
ferent atomic numbers A. As already discussed there, 
this could represent an important difference compared to 
PSB because an increasing threshold energy may allow 
the nucleus to propagate over longer distances. For the 
cross sections crA.i{e') above threshold we used the same 
parametrization as PSB. 



II. ENERGY LOSS RATES 

Heavy nuclei are attenuated basically by two processes: 
photodisintegration on the diffuse photon backgrounds 
and creation of pairs (l6[ ^ . For energies above 
10^*^ eV, it is the CMB which mostly contributes to the 
photodisintegration process, whereas at lower energies 
the infra-red background provides the main source of 
opacity. 

Pair production occurs at a threshold energy of 2me 
for the photon in the rest frame of the nucleus, and gives 
an important contribution only for the interaction with 
the CMB. We have tabulated the pair production energy 
loss rates from Chodorowski et al. j2|] and treat them as 
continuous losses pof . 

The rate for photodisintegration is given by po| : 
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where A is the atomic mass of the nucleus and i is the 



We include the contributions from three different com- 
ponents of the photon background: the first is given by 
the infra-red photons emitted by galaxies and extends 
from ~ 3.0 x 10"^ eV to ~ 0.33 eV. We used the new es- 
timates obtained from the emissivity of the IRAS galax- 
ies [|4|. The second one is the CMB, extending from 
~ 2.0 X 10-6 eV to ~ 4.0 x lO'^eV, and the third is 
the universal radio background (URB) extending from 
~ 3.0 X 10-9 eV to 3.0 x 10"^ eV. Due to Galactic con- 
tamination the latter can not be measured directly below 
~ 1 MHz, however, we verified that even the highest the- 
oretical estimates from summing over the contributions 
from normal and radio-galaxies [ p?! result in a negligi- 
ble contribution to photodisintegration at the energies 
of interest (< lO^^eV). Finally, we neglected the optical 
background because, as can be seen in Fig. 1 in Ref. [ p2| , 
it has no significant effect. 

In a photodisintegration event the changes in energy, 
AE, and atomic number, AA, are related by AE/E = 
A A/ A. Thus, the energy loss time due to photodisinte- 
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gration is given by A/Rcs,a, where 
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In Fig. l] we show the energy loss time due to single- 
nuclcon, double-nuclcon and multi-nucleon emissions in 
the combined CMB, infra-red and radio background, for 
different atomic numbers. We note that at energies above 
10^*^ eV the heaviest nuclei start to disintegrate more 
quickly. In addition, at these energies, the multi-nucleon 
emission becomes more important compared to one or 
two nucleon emission. Note also that the energy losses 
for ''He (shown as a solid line in Fig. |l]) do not include 
photopion production, that become significant for ener- 
gies > 1.5 X 10^0 eV. 



numerical resolution limits us to Ic ^ 0.008L. We use 
Ic — 0.01 Mpc. The magnetic field modes are dialed on a 
grid in momentum space according to this spectrum with 
random phases and are then Fourier transformed onto the 
corresponding grid in location space. The r.m.s. strength 
B is given by = dkk"^ (_B^(fc)). The simulations 
have been performed for two different strengths of the 
magnetic field: a weak field corresponding to 10~^^ G 
and a strong field corresponding to 2 x 10~^G. 

We injected 6 • 10^ Iron nuclei at the source. The equa- 
tions of motion in the presence of the magnetic force and 
the continuous energy loss due to pair production are 
solved and at least every 0.01 Mpc the nucleus is tested 
against photodisintegration and photopion production, 
by using the rates determined as described in the previ- 
ous section. 

We keep track of each individual secondary nucleus 
and each time such a particle crosses one of the spheres 
of a given radius around the source, arrival direction and 
energy are registered as one event on this sphere. Energy 
loss processes and deflection are treated equally for all 
produced secondary nuclei and nucleons. In the diffusive 
regime each trajectory is followed for a maximal time of 
10 Gyr and is abandoned if the particle reaches a linear 
distance from the source that is twice the distance to the 
furthest sphere. 



IV. RESULTS 
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FIG. 1: The energy loss time vs energy for photodisintegra- 
tion on the combined CMB, infra-red and radio background. 
The solid line is for Helium nuclei, the dotted line for Carbon, 
the dashed line for Silicon and the dot-dashed line for Iron. 



III. NUMERICAL SIMULATIONS 

We use the same numerical approach used in earlier 
publications ^ |2^, putting a single source at the 
center, and we register all nuclei and nucleons arriving 
on 10 "detector" spheres surrounding the source at radii 
scaled logarithmically between 1.5 and 50 Mpc. If not 
indicated otherwise, we also assume that the source in- 
jects Iron nuclei with a spectrum extending up to 
~ 10^2 eV. 

We assume a homogeneous random turbulent mag- 
netic field with power spectrum {Bik)"^) cx fc"^ for 
27r/i < k < 2tt/Ic and {B^{k)) = otherwise. We use 
TLB = —11/3, corresponding to Kolmogorov turbulence, 
in which case i, the largest eddy size, characterizes the 
coherence length of the magnetic field. For the latter we 
use L ~ 1 Mpc, corresponding to about one turn- around 
in a Hubble time. Physically one expects Ic ^ i, but 



We started our simulations injecting a distribution of 
Iron nuclei following an E~'^ power law up to 10^^ eV. We 
then followed their disintegration history and kept track 
of all secondary nuclei produced. We were thus able to 
evaluate the chemical composition of detected nuclei at 
any given distance. 

In Figs. ^ and ^ we show the chemical composition of 
particles detected at three different distances for a mag- 
netic field of 10~^2 Q and 2 x 10"^ G, respectively. Re- 
sults are expressed as integral energy spectra n{> Eth) of 
nuclei of mass A detected above Eth, as a function of A. 
We normalized this quantity to the number nFc(> ^-th) 
of (Iron) nuclei emitted above the same energy. 

In both figures one can see the increasing fraction of 
light nuclei with increasing distance, and the progres- 
sive disintegration of Iron nuclei. In the case of strong 
fields shown in Fig. ||, heavy nuclei are considerably de- 
flected, which implies that the propagated pathlength be- 
fore reaching a given linear distance d from the source is 
much larger than d, the difference being more important 
for high Z and low E nuclei (see spectra below). Diffu- 
sion increases the number density of these nuclei due to 
their increased local residence time, but it also increases 
their probability of photodisintegration at a given dis- 
tance. Therefore a strong magnetic field, i.e. such that 
some UHECR enter a diffusion regime, not only modifies 
the energy spectrum, it also modifies the chemical com- 
position at a given distance, with respect to the case of 
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FIG. 2: Number of nuclei relative to Iron detected above 
4 X 10^^ eV at three different distances from the source as 
a function of their atomic mass for a field B — 10~^^ G. 
Solid, dotted and dashed curves correspond to distances 
d = 1.5 Mpc, 7.1 Mpc, and 50 Mpc, respectively. 



FIG. 4: Relative chemical composition above 10^^ eV as a 
function of the distance for B — 10~^^ G. 




FIG. 3: Same as in Fig. Q but for B = 2 x 10"* G. 



rectilinear propagation (small deflection limit in a weak 
magnetic field) . Further effects of the interplay between 
magnetic diffusion and energy losses will be shown below. 

To show the photodisintegration histories for different 
nuclei, we plot in Figs. ^ to ^ the relative abundances 
fi{d) of various atomic species as a function of distance 
d for different threshold energies, where fi{d) is defined 
as 



ni{d) 



(3) 



and Hi (r) is the number of nuclei of species i detected at 
distance d from the source. 

As expected. Iron dominates the chemical composi- 
tion at small distances, whereas only protons are left for 
very large distances, in agreement with previous stud- 
ies on Iron nuclei propagation (in the absence of a mag- 
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FIG. 5: Same as Fig. || but for B = 2 x 10"* G and £ > 
10^^ eV. 



netic field). However Figs. ^ and ^ show that UHECR 
above 10^" eV cannot be predominantly Iron at distances 
larger than ~ 10 Mpc |^ when propagating in a field of 
strength ~ 2 • 10"^ G. In the case of a weak magnetic 
field, this component can survive with a fraction ^ 10% 
at all energies at distances up to ~ 50 Mpc. Again the 
effect of the magnetic field is due to diffusion which in- 
creases significantly the effective propagated distance for 
a given linear distance. 

The effect of diffusion also becomes apparent by com- 
paring the relative abundances of heavy nuclei in Figs. ^ 
and |5|: The case of a strong field shows an enhancement of 
the relative abundances of Si, O and C by up to a factor 
100 with respect to the weak field case. At higher thresh- 
old energy this effect becomes negligible because nuclei 
are no longer in the diffusive regime. Light elements such 
as Helium are continuously produced by photodisintegra- 
tion of heavier nuclei, they reach a maximum relative 
abundance, which we found to be around 1%, then they 
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FIG. 6: Same as Fig. for a weak magnetic field but for FIG. 8: Same as Fig. H for a weak magnetic field but for 
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FIG. 7: Same as Fig. g for a strong magnetic field but for 
£ > 4 X 10^^ eV. A fine for a given species that stops at 
distance < 50 Mpc means tfiat there is no particle left in tfiat 
species at greater distances in the simulation. 



FIG. 9: Same as Fig. |5| for a strong magnetic field, but for 
E > 10^" eV. A line for a given species that stops at distance 
< 50 Mpc means that there is no particle left in that species 
at greater distances in the simulation. 



quickly disappear, reducing their abundance to 0.01 or 
0.001% at a distance of 20 Mpc. 

One should note that the above figures are sensitive to 
the initial maximum injection energy. In effect, here this 
energy E'max = 10^^ eV, which means that one cannot 
detect (secondary) protons with energy E > i?max/56 ~ 
1.8 • 10^° eV. If the maximum injection energy is low- 
ered, say i?max — 10^^ eV, then one would not see pro- 
tons with energy i? > 1.8 • 10^^ eV, and consequently, in 
Figs. the chemical composition would be domi- 

nated by Iron nuclei at all energies. In Figs. |^J^ with 
threshold E'th = 10^^ eV, the proton domination would 
be reduced. Furthermore, in Figsj^J^, the composition 
would becomes dominated by intermediate mass nuclei 
at distances > 15 Mpc. 

These figures also depend on the energy spectrum in- 
dex chosen. Indeed, it is easy to see that if the energy 



spectrum of primary nuclei A'^p(> E) oc i?^^", then if 
all nuclei above energy E are photodisintegrated in sec- 
ondary protons of energy EjA (and above), the number 
ratio of secondaries A^s(> EijA) to primaries at the same 
energy iVp(> EjA) reads 7Vs/-/Vp = A^"" (all other losses 
neglected). Therefore, depending on the spectral index 
a (taken here as a = 2), the secondary flux has more or 
less importance compared to the primary flux. For hard 
spectra a < 2, the secondaries tend to dominate, while 
the reverse is true for a > 2. 

This latter statement is obviously modified in the pres- 
ence of a strong magnetic field, since particles of a same 
energy but different mass have a difference magnetic 
rigidity. As a consequence, at a same energy, high Z 
particles (in our case, primaries) may be diffusing and 
their local density increased while low Z particles (e.g., 
here secondary protons) may be non-diffusing and their 
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local density not increased. In a strong magnetic field, 
for a hard injection spectrum a < 2, one may thus see dif- 
ferent regimes, in which either the secondaries dominate 
(low energy, where both protons and Iron nuclei diffuse, 
or high energy, where both protons and Iron nuclei do 
not diffuse) , or the primaries dominate (when protons do 
not diffuse but Iron nuclei of the same energy diffuse) . If 
a > 2, then secondaries give a subdominant contribution 
in all cases. 

To further investigate the diffusion problem and pho- 
todisintegration processes we studied the energy depen- 
dence of the average mass and the observed spectra at 
different distances from the source. Figs. |l^ and |ll| show 
the average detected logarithmic nucleus mass log A, as 
a function of energy, for two different distances from the 
source. The sudden change of the plots at an energy 
around 2 x 10^° eV is also due to the maximum injec- 
tion energy which translates here for a maximum proton 
energy ~ 1.8 x 10-^° eV. 

In these figures, we see that at low energies < 10^" eV 
the average composition is more strongly dominated by 
Iron nuclei in the strong magnetic field case than in the 
weak magnetic field case. This is an effect of diffusion, 
as before, which increases the local density of diffusing 
particles vs that of non-diffusing particles. While Iron 
nuclei of energy ^ 10^° eV diffuse ini3~2-10^*G, pro- 
tons of the same energy do not diffuse, hence the effec- 
tive enhancement of Iron nuclei with respect to secondary 
protons. Here as well, note that the above conclusion de- 
pends on the spectral index chosen. If the spectrum is 
hard {a < 2) then the importance of the secondary pro- 
ton is increased with respect to that of the primary nuclei 
flux, and the above effect is reduced. 

At higher energies, an opposite effect happens, i.e., the 
composition is lighter for a stronger field, because pho- 
todisintegration is more important than at low energies 
and increases with the larger propagated pathlength in 
stronger fields. For larger distances (see Fig. [Tl] ) a contin- 
uous increase of the average logarithm of A is seen above 
the high energy proton cut-off for both field strengths. 
This is due to the fact that the maximum energy that 
one can detect for a species of mass A increases with A, 
as discussed above, which implies that moving toward 
higher energies we select heavier nuclei. 

Note that here as well these figures depend rather 
strongly on the initial maximum injection energy. More- 
over, the rather large error bars on the average logarith- 
mic mass are not really representative of a gaussian stan- 
dard deviation, since the mass distribution is strongly 
peaked on Iron nuclei and protons. The large error bar 
simply reflects the large mass difference these two peaks. 
This two-peak behavior can be seen in Figs. ||,|| which 
show the distribution in mass of the composition above 
•10^^ eV for various distances: one clearly see in these 
figures that most recorded particles are either protons or 
iron nuclei. This effect is reduced in the case of a strong 
magnetic field (as photodisintegration losses are more se- 
vere due to increased effective length traveled), as shown 
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FIG. 10: Average logarithmic nucleus mass A as a function 
of energy for B = 10"^^ G (diamonds) and B = 2 x 10"* G 
(squares) at a distance d = 1.5 Mpc. 
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FIG. 11: Same as Fig. |T| but at for d = 7.1 Mpc. 



in these figures and by the reduced size of error bars on 
the average log A in Fig. [l0|,[ll| . 

We finally show in Figs. |l2|and[T^ the expected spectra 
at different distances for the two field strengths. Fig. |l^ 
shows the weak field case and can be compared to Fig. 3 
of Rcf. ||2^ (although these authors chose to use a spec- 
tral index a — 3). Figure 12 shows a cut-off at energy 
E ~ 1.5 — 2 ■ 10^" eV that is increasingly pronounced with 
distance, in agreement with previous works |2^, p3[ . 
Figure |l^ shows a characteristic spectral slope due to 
diffusion of nuclei in the magnetic field for energies be- 
low the cut-off, and an almost fiat component at highest 
energies (recovery of the injection spectrum in absence 
of losses, for rectilinear propagation). The fact that the 
transition energy between diffusive and rectilinear propa- 
gation occurs around the cut-off energy ^ 10^° eV is due 
to the choice of the magnetic field strength B ~ 20 nG. 
If the magnetic field were substantially stronger, the in- 
creased length traveled for particles above the cut-off 
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would result in a more pronounced cut-ofF for a same 
distance. 

In Fig. |l^ one also notes the presence of a low energy 
cut-off around E ^ 1.5 ■ 10^^ eV. This is due to the fact 
that the Larmor radius of Iron nuclei at energies around 
10^^ eV and in a magnetic field around 2 x lO^^G is 
about 2 • lO^pc, comparable with Ic, which represents 
the resolution of our numerical simulation. Furthermore, 
we follow nuclei up to a maximum propagation time of 
the order of the age of the Universe, Tmax = 10^° Gyrs; 
particles with high Z and low energies can have a prop- 
agation time larger than Tmax and never reach distant 
shells. This effect represents an additional contribution 
to the low-energy cut-off of spectra on distant shells. 
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FIG. 12: All-particle spectrum observed at distances d — 
1.5,2.3,3.2,4.8,7.1,10.5, 15.5, 23., 33.9, 50 Mpc from right to 
left. The dotted line is for d = 50 Mpc, and B = 10"^^ G. 
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FIG. 13: Same as Fig. |l| but for B = 2 x 10"** G. 



V. APPLICATIONS 

Let us now apply the results of our numerical simula- 
tions to test some models recently proposed to explain 
the origin of cosmic rays. 



A. Iron nuclei from nearby starburst galaxies 

In a recent work Anchordoqui et al. pO| have put 
forward the possibility that cosmic rays above the an- 
kle are essentially heavy nuclei which originate in two 
nearby {d ^ 3 Mpc) sources, the starburst galaxies M82 
and NGC 253, and propagate in a i? ~ 15 nG extra- 
galactic magnetic field which isotropize the arrival direc- 
tions on Earth. They based their analysis on analytical 
estimates of the diffusion coefficient and approximations 
to the photodisintegration losses and angular deflections. 
Our numerical simulations are well suited to improve the 
discussion of their hypothesis, thanks to a more accu- 
rate treatment of photodisintegration processes and to 
a treatment of deflection without approximations. One 
should first note that Fig. ^ shows that starting with a 
distribution of Iron nuclei at a linear distance d ~ 3 Mpc 
with B ~ 20 nG, the average nucleus mass is still high: 
log^ ~ 3 - 4 at ~ lO^^'eV. This suggests that the 
heavy component can survive across this distance; this is 
in agreement with the results of Ref . |^ . 



100 




E [EeV] 



FIG. 14: The average angular deflexion vs. energy for a source 
at a distance d = 3.2 Mpc. 

In Fig. |l^ we show the angular deflection, defined as 
the angle between the source direction and the momen- 
tum of the particle when it is recorded, as a function of 
energy. Here as well, the sudden increase of error bars 
around 2 x 10^" eV is due to the presence of secondary 
protons in the signal; protons with E ~ 10^° eV suffer a 
similar deflection than iron nuclei of energy ~ 2 x 10^"'^ eV, 
which is of order a few degrees. The same angular deflec- 
tion when plotted vs magnetic rigidity R = E/Z shows 
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a much more regular behavior, similar to that shown in 
Fig. [ij up to the size of the error bars. 

This figure shows furthermore that for energies below 
the transition energy, the arrival directions have been 
isotropized as ~ 90° ± 90°. In the high energy regime, 
one recovers a power law behavior oc but we find 
an average angular deflection that is overall a factor ~ 4 
from that given analytically from a random walk argu- 
ment by Waxman & Miralda-Escude [Q, and used by 
Anchordoqui et al. . We believe this difference is due 
to order of unity factors entering the random walk for- 
mula and our convention of defining the coherence length, 
and the range of applicability of the random walk for- 
mula. 

Of interest is the prediction of an anisotropy that 
should be seen at the highest energies, i? ~ 2 — 3 ■ 10^° eV 
as suggested by Fig. We also note that the high- 
est energy Fly's Eye event Fly's Eye event of energy 
i? = 3.2 ± 0.9 X 10^° eV arrived from a direction that 
is ~ 37° away from M82 (see Anchordoqui et al. 
and references therein). By comparing with Fig. this 
event appears to be only in marginal agreement with 
our simulation, but the rather large deflection could be 
explained by a slight overestimate of the energy. Note 
also that the Fly's Eye event is located 98° away from 
NGC 253. 

Two other very high energy events have been reported 
by the AGASA experiment, one with E ~ 2.1 ± 0.6 ■ 
10^0 eV with arrival direction {a, 5) = (19°, -^21°) § 
(equatorial coordinates), the other with E ^ 2> x 10^*^ eV 
and arrival direction (a, 5) = (359°, 22°) (this latter is 
preliminary, see Ref. Q). These two events are located 
at 82° from M82 and 47° from NGC 253 for the former, 
and 86° from M82 and 49° from NGC 253 for the latter. 
For these two events as well the agreement with Fig. |lj 
is marginal, although slightly better respect to the Fly's 
Eye event. 
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FIG. 15: Observed spectrum at 3.2 Mpc for B=2 x 10~**G, 
compared with AGASA data. Injection spectrum oc E~^'^ 
and normalization was fit to the data. 



Finally, in Fig. |l5| we compare the spectrum observed 
at distance d — 3.2 Mpc from a single source for a mag- 
netic field B — 20 nC with the observed AGASA spec- 
trum. It turns out that in order to fit the AGASA spec- 
trum an injection spectrum oc E~^'^ is required. This 
is relatively hard compared to the E^^ injection spec- 
trum usually expected for shock acceleration [Q. The 
hardness of the spectrum required is likely due to the 
interplay between energy losses, existence of secondaries 
and diffusion at low energies. 

Note that it is not realistic a priori to expect that a 
source such as a starburst galaxy would accelerate only 
iron nuclei, and not lighter nuclei. In particular the cos- 
mic abundance of iron would suggest that protons should 
be much more abundant at the same rigidity. Assuming 
that the accelerated spectrum for species i as a function 
of rigidity R can be written drii/dR — Ni{R/ Rq)~°' , 
one finds that the ratio of fluxes of two species i, j at 
a given energy reads FjFj = {Zi/Zj^-^N./Nj. If i 
corresponds to protons, j to iron nuclei, and a — 1.6, 
then Fp/Fpc ~ O.lNp/Npe. If Np/Npc is corresponds to 
the cosmic abundance of iron, then indeed one cannot 
consider iron as the dominant species. However in the 
present scenario, it is assumed that acceleration takes 
place in two steps, first in supernovae shock waves up 
to 10^^ eV then reaccelerated in the galactic wind up to 
10^" eV. It is not clear in this case to what Np/Npc refers, 
but if in a first approximation one considers that it is the 
proton to iron ratio at ~ 10"'^^ eV, this latter is found to 
be of order a few for galactic cosmic rays. This is mainly 
due to the fact that the spectrum of heavier nuclei cos- 
mic rays is generically harder than that of lighter nuclei. 
In that case, indeed the contribution of protons, and for 
that matter, or intermediate mass nuclei, to the energy 
spectrum at injection can probably be neglected in a first 
approximation. 

Finally, one should note that starburst galaxies are ac- 
tive for a finite amount of time: only ^ 10* yrs. Here for 
S ~ 20nG the time delay at £; - 10^^ eV is of order of 
a few 10* yrs. If the high energy part of the spectrum 
has been recorded (in part) by the present experiment, 
then the flux at the lower end of the spectrum should 
be depleted in heavy nuclei, as most of these particles 
would not have had enough time to reach us. It is diffi- 
cult to quantify this effect at present, but it constitutes a 
potential signature of this scenario for future detectors. 



B. Helium nuclei from M87 

In a different scenario, proposed by Ahn et al. |Q, M87 
in the Virgo cluster (located at a distance d « 20 Mpc 
from the Milky Way), is assumed to be the local source 
of UHECR. Indeed the authors showed that one can trace 
back to M87 the 13 events observed above 10^° eV §] if 
the Galactic magnetic field has the structure of a Parker 
spiral and extends to ~ 1 — 2 Mpc. They showed that 
provided the two highest energy events are Helium nuclei 
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and the others protons, all 13 events point back to within 
20° of M87. The importance of the specific magnetic field 
chosen to reach this conclusion was stressed in a note by 
Billoir and Letessier-Selvon ||3^ . 

To see if such a composition is possible we studied the 
relative abundance of Helium as a function of the dis- 
tance from the source. In Fig. |l^ we show the chemical 
composition as a function of distance, assuming that only 
He nuclei are injected in M87 and that the extra-galactic 
magnetic field B = IQ-^^ G. 



Nucleons; 
Helium 



Distance(Mpc) 

FIG. 16: Relative chemical composition as a function of the 
distance for a simulation with a source injecting Helium nu- 
clei, for B = 10-^2 G. 

This shows that even for very weak fields and thus neg- 
ligible deflection, the abundance of Helium nuclei with 
energy above 10^° eV is a factor 100 smaller than the 
nucleon abundance at distances d ~ 20Mpc. The prob- 
ability of observing two helium nuclei out of 11 protons 
in this energy range is thus extremely small. If all the 
events are protons, the convergence in the direction of 
M87 is poor, which makes the model much less attrac- 
tive. To start with a distribution of Iron nuclei would 
make things worse, as can be seen in Fig. |6[ leading to 
Helium nuclei abundances at least 10'^ times smaller than 
the nucleon abundance. 

Although we considered a weak stochastic magnetic 
field while Ahn et al. considered a strong coherent mag- 
netic field structured as a Parker spiral, this should not 
make a big difference in our conclusions. As a matter 
of fact, if one increases the magnetic field strength, the 
effective length traveled for He nuclei is increased, hence 
photodisintegration should be more severe. This scenario 



thus appears fine-tuned in so far as the chemical compo- 
sition is concerned. 



VI. CONCLUSIONS 

In this paper we studied the propagation of a distri- 
bution of heavy nuclei in a stochastic magnetic field, in- 
cluding all relevant energy loss processes. For the propa- 
gation in the magnetic field we used the same numerical 
approach as in Ref. |2^, This approach was here 
generalized to heavy nuclei and their photodisintegration 
processes. One main conclusion of this paper is that a 
strong magnetic field, i.e. such that some UHECRs expe- 
rience a diffusive propagation regime, can strongly mod- 
ify the chemical composition and energy spectrum at a 
given energy with respect to what would be seen in the 
absence of a magnetic field. Rather generically, an in- 
creased magnetic field implies a larger effective length 
travelled, hence a larger photodisintegration probability, 
hence a chemical composition shifted to ligther species. 
As we have argued, the extent of this effect also de- 
pends on the injection spectrum spectral index, and on 
the maximal injection energy. If the injection spectrum 
dn/clE oc then if a > 2 the secondary protons pro- 

duced in photodisintegration interactions do not give a 
dominant contribution in the low energy observed flux. 
The converse is not generally true in the case of a strong 
magnetic field, as the injection spectrum is softened by 
diffusion. 

We applied our results to the discussion of two models 
recently proposed to explain the origin of UHECR. Our 
simulations suggest that the model proposed by Anchor- 
doqui et al. in which UHECR are iron nuclei accel- 
erated in nearby starburst galaxies, is in relatively good 
agreement with the data as far as the energy spectrum 
is concerned. However, it requires a relatively hard in- 
jection spectrum (a ~ 1-6), and the three highest energy 
events from AGASA and Fly's Eye are ^ 40° away from 
the galaxies proposed as sources, in marginal agreement 
with the expected deflection. 

We also showed that for an injection spectrum dom- 
inated by Helium nuclei, the relative abundance of He- 
lium compared to nucleons turns out to be smaller than 
0.01 at distances ~ 20Mpc from the source. This implies 
that the scenario of Ahn et al. |Q, which suggests that 
the UHECR originate from M87 and are deflected in a 
powerful Parker spiral Galactic magnetic field, and which 
requires that the two highest energy cosmic rays (out of 
13 above 10^° eV) are He nuclei, is highly fine-tuned. 
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